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Summary 

This  interim  report  summarizes  the  most  relevant  findings  during  the  first  year  of 
research  on  the  effort  titled  Average  Likelihood  Methods  for  Classification  of  CDMA. 
The  goal  of  the  research  is  to  present  the  development  of  several  analytical  form  of  the 
average  likelihood  for  CDMA  using  combinations  of  BPSK  and  QPSK  symbols.  In  a 
previous  report,  the  case  of  CDMA  using  a  BPSK  code  matrix  and  a  BPSK  data  vector 
was  studied  and  published  in  Technical  Report  AFRL-RI-RS-TR-2014-093.  In  this 
report,  the  average  likelihood  for  different  CDMA  transmission  models  will  be 
developed  using  similar  definitions  and  propositions. 

Decision  theoretic  methods  based  on  average  likelihood  are  the  only  methods  that 
guarantee  the  best  classification  in  noise  according  to  the  classical  Detection  Theory. 
Currently  there  is  a  lack  of  average  likelihood  methods  for  CDMA.  This  research  has 
the  purpose  to  establish  a  foundation  for  new  classification  and  estimation  of  CDMA 
signals. 

Keywords:  DS/CDMA  signals,  BPSK,  QPSK,  average  likelihood  function,  classification, 
estimation. 
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1.  Introduction 


The  characterization  of  signals  in  terms  of  the  modulation  parameters  is  a  field  of 
interest  relevant  to  military  and  software  defined  radio.  Early  modulation 
classification  algorithms  were  based  on  ad-hoc  methods.  These  methods  were 
applied  for  classifying  analog  signals  such  as  AM  and  FM.  The  methods  were 
extended  to  digital  types  such  as  BPSK  and  QPSK.  A  collection  of  these  algorithms 
can  be  found  in  [1].  The  main  benefit  of  these  Ad-Hoc  algorithms  is  in  providing 
classification  rules  that  are  easy  to  implement.  The  mayor  disadvantage  is  their 
inability  to  guarantee  any  degree  of  performance  in  noise. 

But  the  field  of  modulation  classification  underwent  a  significant  change  of  direction 
with  Kim,  Huang  and  Polydoros  papers  [2]  and  [3].  These  papers  applied  classical 
detection  theory  to  MPSK  signals.  Starting  from  a  model,  these  authors  developed  an 
average  likelihood  approximation  for  MPSK  signals.  The  authors  also  demonstrated 
that  some  of  the  heuristic  classification  rules  used  to  classify  signals  were  special 
cases  of  the  Average  Likelihood  Ratio  Test  (ALRT).  With  these  findings,  following 
authors  began  the  development  of  likelihood  methods  for  a  combination  of  QAM  and 
MPSK  methods  as  shown  in  Table  1. 

Table  1  was  presented  in  a  2007  conference  paper  [4].  As  we  can  see,  there  is  a  lack 
of  decision  theoretic  methods  for  CDMA  signals.  In  2012,  the  same  authors 
published  a  paper  on  the  classification  of  CDMA  signals  based  on  its  code  length  [5]. 
The  approach  consists  on  extracting  cyclostationary  features,  i.e.,  peaks  in  the 
spectral  correlation  function,  and  then  training  a  neural  network.  Feature  based 
approaches  started  to  emerge  as  an  alternative  of  Decision  Theoretic  methods  with 
increase  complexity.  As  in  the  case  of  Ad-Hoc  methods,  feature-based  approaches  do 
not  guarantee  performance  and  the  main  weakness  of  this  method  is  the 
dependence  of  the  classification  on  the  training  set.  In  addition,  the  proposed 
method  for  CDMA  signals  can  process  signals  for  a  limited  set  of  code  lengths  due  to 
its  inherent  complexity. 

Developing  decision  theoretic  approaches  for  CDMA  signals  is  a  challenging  task  due 
to  the  large  number  of  unknown  variables  that  require  averaging.  Reference  [6] 
shows  that  it  is  possible  to  express  the  likelihood  in  a  compact  form  and  apply 
further  simplifications.  The  detection  rules  were  applied  the  classification  of  signals 
up  to  code  lengths  of  212. 

During  the  second  year  of  performance  of  this  AFOSR  funded  effort,  the  research 
explored  the  extension  of  average  likelihood  for  CDMA  signals  constructed  from  a 
combination  of  BPSK  and  QPSK  symbols.  The  research  found  that  it  was  possible  to 
develop  average  likelihood  equations  using  previous  findings. 
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The  usage  of  the  average  likelihood  applied  to  CDMA  is  not  limited  to  the  problem  of 
classification.  A  more  ambitious  goal  is  to  be  able  to  estimate  the  spreading  matrix 
and  multi-user  detection.  These  problems  are  considered  for  future  work. 


Table  1  Summary  of  Likelihood-Based  Classifiers 


Authors 

Classifiers 

Modulations 

Unknown 

Parameters 

Channel 

Sills 

ALRT 

BPSK,  QPSK, 

16QAM,  V29, 

3 2 QAM,  64QAM 

Carrier  phase 

AWGN 

Wei  and  Mendel 

ALRT 

16QAM,  V29 

- 

AWGN 

Kim  and  Polydoros 

Quasi-ALRT 

BPSK,  QPSK 

Carrier  phase 

AWGN 

Huang  and 

Polydoros 

Quasi-ALRT 

UW,  BPSK,  QPSK, 
8PSK,  16PSK 

Carrier  phase, 
timing  offset 

AWGN 

Sapiano  and  Martin 

ALRT 

UW,  BPSK,  QPSK, 
8PSK 

" 

AWGN 

Long 

Quasi-ALRT 

16PSK,  16QAM,  V29 

Carrier  phase 

AWGN 

Hong,  Ho 

ALRT 

BPSK,  QPSK 

Symbol  level 

AWGN 

Beidas  and  Weber 

Quasi-ALRT 

32FSK,  64FSK 

Phase  jitter 

AWGN 

Beidas  and  Weber 

Quasi-ALRT 

32FSK,  64FSK 

Phase  jitter,  timing 

AWGN 

Panagiotu 

GLRT.HLRT 

16PSK,  16QAM,  V29 

Carrier  phase 

AWGN 

Chugg 

HLRT 

BPSK,  QPSK,  OQPSK 

Carrier  phase, 
signal  power,  PSD 

AWGN 

Hong  and  Ho 

BPSK,  QPSK 

Signal  Level 

AWGN 

Hong  and  Ho 

HLRT 

BPSK,  QPSK 

Angle  of  Arrival 

AWGN 

Dobre 

HLRT 

BPSK,  QPSK, 

16QAM,  V29, 

3 2 QAM,  64QAM 

Channel  amplitude, 
phase 

Flat  Fading 

Abdi 

ALRT,  quasi-HLRT 

16QAM,  32QAM, 
64QAM 

Channel  amplitude 
and  phase 

Flat  Fading 
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2.  Approach 

The  approach  for  deriving  the  average  likelihood  function  for  CDMA  consisted  on 
expressing  the  average  likelihood  in  a  compact  form.  The  function  can  be  expressed 
as  a  sum  of  product  of  hyperbolic  cosine  functions.  Further  simplification  can  be 
achieved  by  implementing  a  hard  decision  classifier.  In  this  approach,  the  variance 
of  the  probability  of  the  code  matrix  is  reduced  to  a  minimum  value.  This  is 
equivalent  to  classifying  the  CDMA  signal  using  the  most  probable  code  for  a  given 
code  length. 

The  accomplished  work  on  BPSK  codes  and  BPSK  data  considered  the  case  of  chip 
synchronous  signals  with  no  knowledge  on  the  beginning  of  the  sequence.  Both 
balance  and  unbalanced  load  cases  were  considered  under  the  synchronization 
assumption.  A  similar  methodology  will  be  applied  to  general  CDMA  models  using  a 
combination  of  BPSK  and  QPSK  symbols.  The  following  CDMA  transmission  models 
were  investigated. 

•  QPSK  code  -  BPSK  data, 

•  BPSK  code  -  QPSK  data,  and 

•  QPSK  code  -  QPSK  data. 

3.  Development  of  the  Average  Likelihood 
3.1  General  Model 

The  goal  is  to  classify  all  types  of  CDMA  signals  by  their  code  length  L  and  number  of 
users  U  in  the  presence  of  Additive  White  Gaussian  Noise  (AWGN).  A  CDMA  signal  x 
of  size  L*  x  1  is  generated  by  modulating  a  spreading  code  matrix  C  of  size  L*  x  U*, 
with  the  information  symbol  vector  b  of  size  U*  x  1.  Both,  matrix  C  and  vector  b  are 
constructed  from  a  combination  of  BPSK  {+ 1}  and  QPSK  {+ 1,  ±i]  symbols. 

For  simplicity,  we  assume  that  the  energy  per  symbol  per  user  E  is  constant,  i.e.,  a 
balanced  load.  The  code  length  and  number  of  active  users  is  assumed  to  be 
unknown  in  the  model;  therefore,  instead  of  using  the  true  dimensions  L*  and  U*, 
the  model  uses  the  hypothetical  dimensions  L  and  U  as  follows: 

x  —  y/E/L  C  b 

(1) 

The  spreading  matrix  is  a  block  matrix  composed  of  column  vectors  c* with  low 
cross-correlation.  A  metric  that  characterizes  the  spreading  matrix  is  known  as  the 
Total  Square  Correlation  (TSC).  A  modified  TSC  is  defined  in  this  paper  by  [7] 
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(2) 


u- 1  u-i 

t(C)  =  ^  ^\(ctii,c,j)\Z 
i= 0  j=0 

j*i 


An  alternative  definition  of  (2)  for  a  code  matrix  with  either  BPSK  or  QPSK  symbols 
can  be  expressed  in  terms  of  the  Frobenius  Norm  with  L  is  the  number  of  rows  of 
CHC  and  U  is  the  rank  of  CHC: 


T (C)  =  \\CH C\\2F  -  L2  ■  U 


(3) 

The  noise  vector  n  has  zero  mean,  and  covariance  matrix  £  =  N0/ 2  ■  I,  where  /  is 
anLxt  identity  matrix.  The  receiver  has  no  knowledge  of  the  spreading  code,  the 
code  length  or  the  number  of  users.  For  a  given  hypothesis  0-C  —  { L ,  U },  the  receiver 
makes  the  assumption  that  the  code  length  is  L  without  necessarily  implying  that 
L  —  L*.  The  channel  is  modeled  by: 


ye  —  x£  +  n 


(4) 


where  ye  represents  the  received  vector  delayed  by  e  number  of  chips. 


The  likelihood  is  conditioned  on  the  hypothesis  0-C  —  { L ,  U},  the  spreading  matrix  C, 
and  the  information  vector  b  .The  delay  xe  is  unknown  to  the  receiver;  therefore,  the 
comparison  is  done  with  respect  to  x  from  a  hypothetical  model.  Under  the 
assumption  of  no  frame  synchronization,  our  conditional  likelihood  is  given  by: 


A(y\K  =  {L,  U},  C,  b,  e)  =  ^nN^L/2  e 


(5) 


3.2  Average  Likelihood  Function 

The  construction  of  our  likelihood  function  requires  averaging  over  the  variables 

— > 

C,  b  and  e.  The  variable  e  will  be  ignored  in  this  development  due  to  its  lack  of 
mathematical  interest.  For  convenience,  we  define  the  energy-to-noise  ratio  per 
chip  y,  the  unnormalized  signal  s,  and  the  correlator  output  r. 


Y 


E 

N0L’ 


C  b,  R  — 


y 


(6) 
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The  conditional  likelihood  function  of  (5)  is  a  function  of  the  correlator  equation 
RhS  +  ShR  and  the  energy  of  the  signal  SHS  as  follows: 


A(f\K,  C,  b ) 


fe -<*’*>  Y 1/2 


ejY/2(RHS+SHR)~Y  SHS 


(7) 

The  development  is  shown  in  Appendix  A.  In  Detection  Theory,  the  constant  term 
(R,  R)  is  usually  omitted;  however,  the  hypothesis  depends  on  the  code  length  and 
the  receive  vector.  Caution  should  be  exercise  before  trying  to  eliminate  terms  that 
depend  on  R. 

It  was  verified  that  a  classification  rule  developed  under  the  assumption  of  a 
uniformly  distributed  probability  of  the  code  does  not  work.  A  key  observation  is 
that  there  is  a  dependency  between  the  elements  of  the  spreading  matrix.  It  is 
unlikely  that  a  random  selection  of  matrix  elements  would  result  in  a  low 
correlation  matrix.  Therefore,  we  must  selectively  choose  the  spreading  matrices 
based  on  some  code-design  metric.  The  proposed  approach  constructs  a  discrete 
probability  distribution  based  on  our  definition  of  the  TSC.  Our  first  choice  is  to 
propose  a  discrete  probability  that  present  similarities  with  a  Gaussian  distribution: 

e-/?-T  2(c)w(c) 

PfCl  — - — — — 

SC'e-Ww(C0’ 


(r fO  rank(C )  =£  U 
=  ll  otherwise 


(8) 

The  function  w(C)  is  a  weight  that  eliminates  inadmissible  spreading  matrices,  i.e., 
matrices  with  linearly  dependent  codes.  For  perfectly  uncorrelated  codes  C* ,  the 
TSC  metric  is  zero  and  the  probability  achieves  a  maximum  value.  Due  to  the 
invariance  of  the  TSC  to  permutation  matrices,  any  code  matrix  Q*that  is  a 
permutation  of  C *  will  have  the  same  probability.  In  the  coding  theory,  Q *  and  C* 
are  referred  equivalent  matrices  [8], 

The  parameter  /?  is  called  the  precision  of  the  distribution  in  Machine  Learning 
Theory.  It  is  defined  as  half  the  inverse  of  the  variance  of  a  Gaussian  distribution  [9]. 
In  the  case  of  a  Gaussian  distribution,  the  precision  reshapes  the  probability  by 
making  it  steeper.  At  the  same  time,  it  also  affects  the  classification  boundaries,  by 
changing  the  steepness,  thus  converting  the  classification  from  a  soft  decision  to  a 
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hard  decision  problem.  Averaging  over  the  code  coefficients  eliminates  the 
dependency  of  the  likelihood  on  the  spreading  matrix  as  follows: 


A(r\H,  b ) 


eJY/2(RH S+ShR)~y  shs 


(9) 


Finally,  the  dependency  of  the  conditional  likelihood  on  b  is  removed  by  averaging 
over  uniformly  distributed  information  symbols  as  given  by: 


A(r\0-C)  = 


/e-(R,R)\1/2 
No. TL) 


q-JyTA  ( rhs+shr)-y  shs 


(10) 

The  probability  P(b )  depends  on  the  domain  of  b  with  P(b )  =  2~u  for  BPSK 
symbols  and  P(b)  =  2~2U  for  QPSK  symbols.  The  order  in  which  the  expectation  is 
taken  does  not  affect  the  final  result.  In  our  study,  the  average  over  the  information 
vector  is  done  first  for  convenience  as  we  will  see  later  in  the  discussion. 

3.3  Key  Definitions  and  Propositions 

The  inclusion  of  complex  variables  increases  the  complexity  of  the  problem  by  the 
number  of  unknown  variables  when  compared  to  the  BPSK  code/BPSK  data  CDMA 
model.  In  order  to  average  the  conditional  likelihood,  we  provide  the  following  key 
definitions  and  propositions. 

r  -.LXU 

Definition:  The  set  S$  is  the  subset  of  real  antipodal  matrices  C  —  [Cij  j  with: 

u-i 

o-i  =  ^  Cij  ^  0 

7=0 

(11) 

L'X.U 

Definition:  S^  t  is  the  subset  of  real  antipodal  matrices  C  —  [cij]  E  S$  with: 


t2(C)  =  t 

(12) 

Definition:  t  is  the  subset  of  pair  of  matrices  (Z+,Z_)  with  Z+  =  {±l}LxU 

andZ_  =  {+l}ixU  with: 
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U-l 

(a+)i  =  ^  0 

7=0 

U-l 

(5_)i  =  ^]  iz-)i,j>0 
7=0 

/(Z++Z_)  (Z+-Z_)\ 

t  v — 2 — +n — 2 — ; =t 

(13) 

Proposition  1:  The  product  of  a  diagonal  matrix  GLxL  —  diag(gLxl )  and  a  vector 
rLxl  commutes  as  follows: 


G  ■  r  —  diagif )  •  g 


(14) 


Proposition  2:  The  product  of  the  block  matrix 

r  \  Gr  G, 

L-g 7  gJ 

with  diagonal  matrices  GR  =  diag(gR )  and  G/  =  diag(gj)  and  a  block  vector  R2Lx  1 
with  entries 


R 


2Lxl 


rR+r, 


commutates  as  follows: 


GR 


diag(r ,) 
—diag  (rfi) 


diag  (rR)l  \gR  +  gl 
diag (r7)J  -  g,. 


(15) 
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Proposition  3:  The  average  of  a  function  f(C  ■  b)  over  C  and  b  where  C  E  {±l}Lxt/ 
and  b  E  {±l}t/xl  is  given  by: 


Z /(c  ■ g) = fic  ■  i) 

C  b  C 

(16) 


Proposition  4:  Let  b  —  bR  +  i  bj  with  h  G  (±1,  +  I  }c/xl;  bR  —  Re{b)  and  bI  —  Im{b). 
Also,  let  b+  E  {±l)(7xl  and  G  {±l}t/xl,  then 

+  b,))f(pT(iR-b,))  = 

bR  bj  b+  b _ 

(17) 


The  previous  proposition  states  that  for  averaging  purposes,  the  product  b+  and  b_ 
acts  as  two  independent  antipodal  vectors. 


Proposition  5:  The  average  of  the  function  eaTQ  over  g,  where  g  —  (±l)Lxl,  is 
given  by: 


L  —  l 

_L  ^  eaTg  =  Y [  cosh(cfj) 
g={±  i}Lxl  i= 0 


(18) 


Proposition  6:  The  TSC  is  invariant  to  permutations  and  sign-inversions  of  rows 
and  columns. 

Proposition  7:  The  limiting  case  as  /?  -»  oo  of  a  rational  exponential  function 
with  t f  >  0,  di  >  0,  fi  >  0  is  given  by: 

Ef=0  die-rt  d0 

lim - 5-  =  — 

"-“ZiLo  Ke-ftl  f« 

(W) 
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The  proofs  of  these  propositions  are  found  in  Appendix  C. 

3.4  Development  of  a  Compact  Average  Likelihood 

For  each  of  the  three  CDMA  transmission  models  in  section,  we  need  to  find  suitable 
models  that  use  a  development  similar  to  the  already  solved  case,  i.e.,  the  BPSK- 
code/BPSK-data  model. 

The  development  requires  computing  the  correlation  term  ( RHS  +  SHR ),  the  energy 

term  SMS  of  the  conditional  likelihood  of  (10)  and  the  covariance  matrix  of  each 
CDMA  model.  The  covariance  matrix  is  uncorrelated  with  a  common  variance  o 2  for 
all  the  variables. 

3.4.1  Revisiting  BPSK-Code/BPSK-Data  Model 

The  model  for  a  BPSK-Code/BPSK-Data  transmission  model  is  given  by: 


Gr  —  diag(gR ) 

gR  e  (±1}LX1 
bR  e  {il}"*1 
Qr  e  {±1}LXU 


(20) 


or  equivalently, 

t 

0 

R  =  S  +  N 


R  = 


rR 

fi 


S  = 


Gr  0 

o  Gr} 


Qr  0 

LO  Qr 


(21) 


The  imaginary  part  of  (21)  does  not  provide  any  useful  information  and  vanishes  in 
the  correlator  equation.  By  using  Proposition  1,  the  correlation  expression  becomes: 
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RtS  —  gR  D  Qr  bR 
D  =  diag(rR ) 

(22) 

Also,  we  compute  the  energy  term,  which  is  given  by: 

STS  =  bTR  QTRQR  bR 

STS  =  1 TZTZ  1 

(23) 

with  Z  =  Qr  diag(b fi).  The  covariance  matrix  for  complex  signals  is  uncorrelated  as 
shown  below: 

E{NNh}  =  ^  ILXL 

(24) 

The  previous  equations  are  substituted  in  (10).  The  average  likelihood  in  its  raw 
form  is  given  by: 

1 

Krm  =  ^P(ffR)^P(Z)^P(bR)e^(^DZ^-rllzIll2 

'  °  '  9r  z  b 

(25) 

with  P(gR )  =  1/2L.  By  splitting  the  matrix  C  as  the  product  of  G  Q,  we  have  slightly 
modified  the  definition  of  P(Q )  or  P(Z).  Averaging  over  gR  —  diag(G )  doubles  the 
number  of  possible  matrices  Z  whenever  a  row  of  the  product  Z  1  is  zero,  i.e., 
multiplying  zero  by  +1  is  not  different  from  multiplying  zero  by  -1.  We  define  P(Z) 
as: 
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e-/?-r2(Z)  2-ZiCount((z  l).==o) 


P(Z)  £ 


Y,Q'  e~P'r2^  2  ^iCount{(z'  ^r- °)  W(Z') 


fr?^  fO  rank(Z )  =£  U 

w(z)  =  tl  0tw« 


(26) 


The  average  over  bR  is  simplified  by  applying  Proposition  3. 


I>&)  g V2 y  (5r  n  <2r  Sr)-/  £r  Qr<3r  bR  _  ej2y  (Sr  D  Z  l)-y  ||z  T 

Pr 


(27) 


The  average  over  gR  is  simplified  by  applying  Proposition  5. 

_  L  — 1 

n  cosh(1yr2y  (D  a)j) 

<?R  i=0 


(28) 


The  summation  over  Z  is  split  in  the  summation  over  all  possible  a  formed  from 
a  —  Z  1  with  Z  G  all  possible  TSC  terms  and  all  matrices  elements  of  the  set 
Sfa.  Matrices  in  the  set  S^t  have  the  same  the  energy  term  and  TSC  according  to 
Proposition  7. 


I-III 


t  ZE.  S-pi 


(29) 


The  resulting  likelihood  is  a  compact  form  given  by: 
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A(r\K)  = 


/  e-(R,R)  \2 


L— 1 

a((3,  a)e~Y^a^2  J” |  cosh(A/27  (D  a)j) 

a  i=0 


a  —  Z  1 


where 


(30) 


a(/U)  ^  P(Z). 


t  ZESt, 


(31) 


The  coefficient  a(/?,  a )  in  (31)  implies  that  detecting  CDMA  signals  using  the 
provided  stochastic  model  requires  a  characterization  of  all  possible  code  matrices 
for  a  given  code  length.  This  fact  can  be  seen  as  a  weakness  of  the  average  likelihood 
approach,  since  the  initial  assumption  does  not  include  the  knowledge  on  the  code 
matrix.  In  practice,  the  user  would  need  to  have  partial  knowledge  on  possible  codes 
used  before  using  equation  (30). 

3.4.2  BPSK-Code/QPSK-Data  Model 

The  model  for  BPSK-Code/QPSK-Data  is  given  by: 


rR 

r,. 


\Gr 

0 

Qr 

0 

bR 

_L 

nR 

0 

Gr. 

0 

Qr. 

■bi. 

1 

n,. 

Gr  —  diag(gR ) 

gR  e  (±1}LX1 

bR  +  5  b[  E  (±1,  +S)Uxl 
Qr  e  { ±1}LXU 


(32) 


The  model  is  impractical  for  developing  a  compact  average  likelihood  function. 
Instead,  an  equivalent  and  more  convenient  model  is  used  for  calculating  the 
average  likelihood: 
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rR+r, 

s  = 

\Gr  0  1 

oi 

bR  +bI 

N  = 

nR  +  n7 

l-(rR-r7)J 

Lo  gr\ 

Lo  qr\ 

-(bR-bj). 

>  1  ¥ 

~(nR  -  n;). 

R  =  S  +  N. 


(33) 


The  equivalent  model  is  a  linear  combination  of  (32).  In  other  to  simplify  the 
development,  we  notice  from  Proposition  4  that  the  summation  over  the  block 
vector: 


bR  +  b, 


(bR  -  b,)\ 


(34) 


acts  as  an  independent  variable  vector  with  binary  antipodal  symbols.  Also,  the 
covariance  matrix  for  complex  signals  is  given  by: 


E{(nR  +  n,)(nR  +  n ,)H}  =  —  ILxL 


(35) 


Only  half  of  the  equation  variables  in  (33)  provide  the  necessary  information  for 
computing  the  average  likelihood.  It  is  possible  to  reuse  average  likelihood  for 
BPSK-code/QPSK-data  CDMA  in  (30)  and  (31)  by  simply  replacing  D  with: 

D  —  diag(rR+rI') 

(36) 


The  BPSK-code/QPSK-data  CDMA  model  splits  the  transmitted  signal  in  two 
orthogonal  components.  The  average  likelihood  recombines  the  real  and  imaginary 
components  of  r  in  (32)  in  a  single  real  vector  rR+r,  to  form  a  model  similar  to  the 
BPSK-code/BPSK-data  model. 
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3.4.3  QPSK-Code/BPSK-Data  Model 

Similar  steps  can  be  used  to  derive  the  formulas  for  all  CDMA  cases  using  BPSK  and 
QPSK  symbols.  The  model  for  QPSK-Code/BPSK-Data  is  given  by: 


rR 

Lr;J 


r  Gr 

- Gi 

Qr 

Qi 

bR 

1 

nR 

[g, 

Gr\ 

[ ~Q , 

Qr. 

.0. 

1 

■W/. 

Gr  —  diag(gR ) 

G,  =  diag{g ,) 

gR  +  igIe{±l,±i}Lxl 

bR  G  (±l)ixl 

Qr  +  IQ,e{±1,±i}LxU 


(37) 


This  model  is  also  impractical  for  developing  a  compact  likelihood  function.  Instead, 
the  equivalent  model  from  Appendix  B  should  be  used: 


rR  +  r, 


Gr  G, 

Qr  ~  Qi 

—Gi  Gr 

+  Qi)  \ 

N  = 


nR  +n, 

~(nR  -«/)■ 


R  =  S  +  N 


(38) 


The  correlator  term  R1  S  can  be  expressed  in  terms  of  two  independent 
variables  g+  —  gR  +  gi  —  {+l)ixl  and  g_  —  gR  —  gi  —  {±l)ixl.  These  vectors  will 
be  averaged  using  Proposition  3. 


RtS  = 


9+ 

T 

-9- 

diag(j{) 


diag(rR ) 
—diag(r,)_ 


Qr  Qi 
Qr  +  Qi. 


(39) 


The  orthogonal  relation  in  (40)  is  derived  from  the  fact  that  GRG,  =  0  and 
GrGr  +  G]  G,  =  I. 
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—  ^2Lx2L 


(40) 


The  proof  is  shown  in  the  Appendix  D.  Thus,  the  energy  of  the  equivalent  model  is 
given  by: 


STS 


Qr  Qi 

h 

l-(Q*  +  Qi)\ 

dr 

2 


STS  = 


(41) 


where  Z+  =  —  (QR  +  Q ,)  ■  diag(bR )  and  Z_  =  ( QR  —  Q ,)  ■  diag(bR )  with  both  Z+ 
and  Z_  G  (±l)Lxt/.  These  matrices  can  be  treated  independently  as  long  as  Z+  =£  Z_. 
For  case  where  Z+  =  Z_,  the  code  matrix  Q  is  real  and  the  model  becomes 
equivalent  to  the  BPSK-code/BPSK-data.  It  is  possible  to  reuse  average  likelihood 
for  BPSK-code/QPSK-data  CDMA  in  (30)  and  (31)  by  simply  replacing  D  and  Z  with: 

_  diag(r diag(rR )  ' 

Ydiagirn)  —diagifj). 


Z  = 


F-l 

z_ 


(42) 


Proposition  4  can  be  applied  to  the  block  vector  [g+,g-]T  in  (39).  Note  that  the 
redefining  Z  in  (42)  does  not  affect  the  parameter  P(Z).  The  new  matrix  ZTZ  has  L 
rows.  The  Frobenius  Norm  contains  an  additional  term:  x 2(Q)  =  t2(Z)  + 

+l|z+z_  —  zZz+\\p. 
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)-/?-(t2(z)+||z: 


u-zlz,  t)w(z) 


Pqpsk(Z)  —  /  I,  T  n2\ 


w( 


[0  rank(Z)  ^  U 
rank(Z )  =  U 


The  resulting  likelihood  is  a  compact  form  given  by: 

1 

~(R,R)\  2  , 


r-ylia-IF 


tt-p  (2_ 


L—l 


(43) 


J- [  cosh(72y  ((r/)i(a+)t  -  (4)i(a-)i))  cosh(72y  ((r*)i(a+)t  -  (r/)t(a-)f)) 


t  =  0 


a+  =  Z+  1 
a_  =  Z_  1 


(44) 


where 


a(/U+,a_)  =  £  ^  P(Z) 

t  (z+,z_)es3+a_  t 


(45) 


3.4.4  QPSK-Code/QPSK-Data  Model 

Finally,  we  can  derive  the  average  likelihood  for  the  CDMA  transmission  model 
using  QPSK  and  QPSK  symbols.  The  model  is  given  by: 
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rR 

fi- 


\Gr 

- Cl 1 

Qr 

Qi 

bR 

nR 

[g, 

Gr 

l-Qi 

Qr. 

b,. 

1 

-n,. 

Gr  =  diag(gR ) 

G,  =  diag(g ,) 

gR  +  igi  e  {±i,±t}Lxl 

bR  +  I  bj  G  {±1,  +i)yxl 
Qr  +  IQ,E{±  l,±i}LxU 


(46) 


The  equivalent  model  provides  a  more  convenient  way  of  developing  the  average 
likelihood  equation. 


rR  +  r, 

,  s  = 

Gr  Gj 

Qr 

Qi 

b+~ 

N  = 

nR  +  n7 

-G,  Gr_ 

-Qi 

Qr. 

-b__ 

>  1  * 

_-(nR  -  ni). 

R  =  S  +  N 


b+ 

b- 


bR  +  b, 
bR  —  K/. 


(47) 


The  correlator  equation  can  be  expressed  in  terms  of  new  variables:  g+  —  gR  +  gIt 
g~  =  gR-  9i,  b&  =  - diag(b+ )  b_,  Z+ =  QR  ■  diag(b+ )  +  Ql  ■  diag(bJ),  and 
Z_  —  Qr  ■  diag(b+ )  —  Q{  ■  diag(b_ ). 


RtS 


9+ 

T 

diag(r — diag(rR ) 

I  0 

Z_ 

-9- 

_diag(rR )  diag(r j)  . 

.0  diag(bA)_ 

Z+. 

Then,  the  energy  is  given  by: 


STS  = 


(48) 


(49) 
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with  bothZ+  andZ_  G  { ±l}LxU . 


The  development  of  the  likelihood  is  almost  similar  to  the  QPSK-code/BPSK-data 
except  for  the  introduction  of  a  new  vector  bA.  We  begin  our  development  by 
averaging  over  by  applying  Proposition  3. 


I>T)  (g  D  F  Zl)-y  \\zif  _  eJty{gRDFzT)-Y  \\Z1 

b+ 


diag(f —diag(rRy 
_diag(rR )  diag(r . 

I  0  I 

F  = 

.0  diag(bA)_ 


(50) 


Then,  we  proceed  to  average  over  g  —  [gR,  cjj]  is  simplified  by  applying  Proposition 

5.  "  . . 

2L-1 

P(gR)6'f^Y  (.Sr  d  f  z T)— y  \\z  i||  _  e-y  ||zT||  J- 1  cosh  (^■s/2y  (D  F  Z  1  )  J 

9r  i=0 


(51) 


The  next  step  is  averaging  over  bA.  This  is  given  by: 

_  2L-1 

^  P(pA)e~Y  Hz  ^11  J” |  cosh  (y2y  (D  F  Z  l)  J  = 


1  =  0 


L—l 


^  />&)  e  7  l|z  111  Y\  cosh  (V^T  (C^/)i  (Z+l).  -  (rR)i  (bA)  ,(Z_l) .)) 
SA  «=o 

cosh  (^27  ((rR)i  (£+  l)f  +  (r7)i  (*A)f(Z_T) .)) 


(52) 
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The  product  of  cosine  can  be  expressed  as  a  sum  of  hyperbolic  cosines  as  follows: 


cosh(y2y  ((f>)i  (Z+l).  -  (fR)j  (hA).(Z_l).)j 


cosh (^[2y  ((rR);  (Z+  l).  +  (r7)7  (bA).(Z_l).))  = 

^cosh  ((rR  +  r7)7  (Z+l).  -  (hA).(rR  -  r7)7  (Z_l).) )  + 
^cosh  (^2/  ((rR  -  r7)7  (Z+l).  +  (hA).(rR  +  r7)7  (Z_l).)) . 


(53) 


The  equation  is  further  simplified  by  applying  the  trigonometric  identity  in  (54). 


-cosh(a  +  d  ■  b)  —  cosh(a)  cosh(b). 

d=±l 


(54) 

Finally,  the  average  over  Z+  and  Z_  is  expressed  in  terms  of  the  set  Z  G 
Note  that  the  redefinition  of  Z  in  (46)  does  not  affect  the  probability  P(Z)  for  the 
same  reasoning  that  was  already  discussed  in  the  previous  section.  See  equation 
(43). 


A(R  |4f)  = 


f  Q  (R>R)  \  2  ^—1  ^—1 
v  /  a+  a.- 


,-Y  l|a+llz-y  lla-ll2 


L— 1 


n(  -cosh  (^2 y  ((rR  +  r7)7  (5+)*))  cosh  (^2 y  ((rR  +  r7)7  (a_)7))  + 

i= 0 

^cosh  (^2 7  ((rR  -r7),  (a+)7))cosh  ((rR  +r7)7  (d_)7))j 


(55) 
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3.4.5  Summary  of  Average  Likelihood  Formulas  for  CDMA 
The  average  likelihood  for  balanced  CDMA  signals  is  given  by: 


A(r  m 


L—l 

J- 1  coshC^/^y  (D  d)i) 

i= 0 


a  —  Z  1 


where 


«(/?,  a)  =  ^ 


I  w 


t  ZeS^, 


Table  2  Average  Likelihood  Parameters 


Code/Data 

Matrix  D 

Matrix  Z  G  S^T 

BPSK/BPSK 

D  —  dia,q(fR ) 

Z  =  ZR 

BPSK/QPSK 

D  —  diagifn+Vj) 

Z  =  ZR 

For  QPSK/BPSK  and  QPSK/QPSK  data,  the  likelihood  coefficients  are  given  by: 

I  P(Z) 

t  (Z+,Z-)ES^+^_t 


«(/?,  aR,  a, )  =  ^ 


Z  = 


P+1 

Z_ 


the  average  likelihood  equations  for  QPSK/BPSK: 


(56) 
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A(f  m  = 


a_)  e-yll«+llz-ylla-ll2 


L— 1 

P~[  cosh(y2y  ((r/)i(a+)i  -  (rR);(a-)i))  cosht^Zy  ((rR);(a+)i  -  (r/)i(a-);)) 

i= 0 


a+  =  Z+  1 
a_  =  Z_  1 


(57) 


and  QPSK-code/QPSK-data: 

1 

/  q  ~{R,R)  \  2  i  i 

A(?|Jf)=  osw  2j2j«w.<M2) 


-y  lla+ll2-y  IIS-II2 


L— 1 


]^~[  ^icosh(72y(rfl  +  f/)i(a+)i)  cosh(y2y(rR  -  r/)i(a_)i) 

-  ^cosh(A/2y(rJ?  -  f/)i(a+)i)  cosh(A/2y(rR  +  r,)i(a_)i)) 


a+  =  Z+  1 
a_  =  Z_  1 

(58) 


3.5  Unbalanced  Load  Classifier 

The  development  of  a  Quasi-Average  Likelihood  function  is  similar  for  all  four  CDMA 
models.  It  is  called  quasi-average  because  the  rule  would  require  the  estimation  of 
the  energy  of  each  CDMA  signature  instead  of  an  average.  A  strictly  average 
likelihood  over  all  possible  signal-energies  per  user  would  result  in  a  complex  sum 
of  product  of  error  functions.  The  models  require  substituting: 


7  + A 

0  1 

bR  +  bj 

— >  — > 

diag(bR  +  b,) 

0 

1  +  5 

0 

/  + aJ 

bR  —  b[. 

— ^ 

0 

diag(bR  —  b7). 

-1  +  5- 

(59) 
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where  the  energy  per  user  Et  is  given  by: 


Et  = 


(V^o(l  + 


2 


(60) 


■4 

E0  is  a  reference  energy,  A  =  diag(5)  and  8  is  a  vector  that  contains  the  deviations  in 
amplitude.  In  other  words,  the  parameter  Z  1  is  substituted  by  Z  (1  +  8 )  likelihood 
equation  in  (30)  as  shown  next. 

^  P(Z)e“ylla+?ll  J~ |  cosh  (^f2y  (p{a  +  f)).) 

All  a  All  r  ZeS^T  i= 0 

a  —  Z  1 
t  =  Z8 

(61) 

3.6  Hard  Decision  Model 

A  simplification  of  (30)  and  (31)  was  proposed  in  [10].  The  likelihood  is  computed 
for  the  case  when  the  precision  parameter  approaches  to  infinity.  The  procedure  is 
found  in  Appendix  (B). 


lim  «(/?,  a)  — 

/?->  CO 


2~Count(ai==0 )  I  r>_ 
_ I  Q.^mir 

s4t„ 


if  a  ^  a 
if  a  —  a 


(62) 


lim  a(/?,  a+,a_) 

/?->oo 


_  1  2_Count((®+)i==0)2_C’ount((a_)i==0) 


>tn 


2-Count((^+)i=-0)2_Count((®-)i  =  =  0)  |£_,  _  ^ 


if  (a+,  a_)  ^  (a+,  a*  ) 

i/  (d+,d_)  =  (d;,di)' 
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(63) 


4.  Status  of  the  Research 

During  this  second  year  of  this  effort,  the  research  has  concentrated  in  expanding 
the  theory  to  cases  of  CDMA  signals  using  BPSK  and  QPSK  symbols.  The  next  task 
will  include  simulations  and  the  preparation  of  a  conference  paper.  The  research 
will  move  to  estimation  problems  in  CDMA  using  the  average  likelihood  as  shown  in 
Table  3. 

Table  3  In-House  Schedule 


Tasks 

Development  of  the  Likelihood 
BPSK-BPSK  Scheme 
Balanced  Energy 
Simulations 
Documentation 
Development  of  the  Likelihood 
BPSK-BPSK  Scheme 
Balanced  Energy 
Simulations 
Documentation 
Estimation 

Estimation  of  the  Code 
Multiuser  detection 
Simulations 
Documentation 


FY-13 

1QTR 


4QTR 


3QTR 

4QTR 
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5.  Conclusions 


It  is  a  known  fact  that  average  likelihood  methods  such  as  ALRT  and  MAP  provide 
the  best  performance  in  noise.  This  means  that  any  other  classification  method,  i.e., 
neural  network  or  ad-hoc  cannot  offer  better  performance  in  noise. 

In  the  case  of  CDMA,  computing  the  average  likelihood  function  from  its  complex 
formulation  in  expression  (10)  turns  quickly  into  an  intractable  problem.  However, 
this  report  shows  that  it  is  possible  to  develop  a  simplified  and  compact  the  CDMA 
average  likelihood  for  combinations  of  BPSK  and  QPSK  symbols.  This  form  consists 
of  a  sum  of  likelihood  coefficients  a  and  products  of  hyperbolic  cosine  functions.  All 
the  studied  cases  of  CDMA  transmission  models  using  a  combination  of  BPSK  and 
QPSK  symbols  can  be  expressed  in  form  similar  to  the  previously  studied  BPSK- 
code/QPSK-data. 

Expressing  the  likelihood  of  CDMA  signals  in  a  similar  form  facilitates  the 
implementation  of  the  algorithms  and  its  validation.  No  significant  issues  are 
expected  from  these  simulations,  since  the  BPSK-code/BPSK-data  has  been 
successfully  tested. 
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7.  Appendices 
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A.  Conditional  Likelihood  for  CDMA 


In  order  to  compute  the  conditional  likelihood,  we  define  the  following  parameters 
and  variables: 


Signal  to  noise  ratio  per  _  E 

chip  Yc  ~  N0L’ 

Unnormalized  signal  s  —  C  -b, 


CDMA  transmit  signal 


Input  to  the  correlator  y 

Under  the  assumption  of  AWGN,  the  conditional  likelihood  can  be  expressed  as  a 
Gaussian  multivariate  distribution  of  a  continuous  vector  y.  Then,  we  proceed  to 
substitute  them  in  the  likelihood  function: 


,  — 1  (y-x)H£  x(y-x) 

A(y\H  =  { L,U},C,b )  = 


A(y\J-C  —  {L,  U},C,b )  = 


\^r-jEjLCb}  r-jEjLCb 

2N0/2 


(nN0)L/2 


w0- 


A{y\K  =  {L,U}X.b)  =  -^A^e  ’ 


A^rHr+  If  J^r(rH  Cb+{Cb]r  )-j{Cb)  Cb 


N0 


A{y\K  —  {L,U},  C,b)  = 


(jiNq)1/2 


-rHr/2  +^[yj2  ( fH  C  b+(c  b)H  f  )-y{Cb)H Cb 


The  final  result  is: 

A(y\H  =  {L,  U},C,b)  = 


__  1/2 

e  {r’r)\  Jyj2(fH  Cb+(Cb)Hr)-y{Cb)H Cb 


(nN0)1 


.  □ 
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B.  Equivalent  Models  for  CDMA 


The  standard  transmission  model  for  any  CDMA  transmission  model 

r  —  G  Q  b  +  n 


can  be  expressed  in  terms  of  real  matrices  as  follows: 


rR 


gr  g, 

Qr 

Qi 

bR 

_l_ 

nR 

-G,  Gr 

Y-Qi 

Qr. 

M. 

i 

.n,_ 

where  <3  =  Qr  +  n  Q,  is  the  code,  b  —  bR  +  S  bl  is  the  data,  G  —  GR  +  I G,  is  a  complex 
diagonal  matrix  that  rotates  the  symbols  in  the  complex  plane,  r  —  rR  +  n  fj  is  the 
correlator  output  at  chip  level  or  the  input  of  to  the  correlator  at  frame  level  and 
n  =  nR  +  I  n;  is  the  noise  signal.  Depending  of  the  CDMA  scheme,  Q  and  b  are 
either  BPSK  (real)  or  QPSK  (complex)  signals.  In  this  report,  we  assume  that  QPSK 
symbols  belong  to  the  set  (+1,  ±n). 


It  was  found  that  a  development  similar  to  the  BPSK  code/BPSK  data  scheme  can  be 
used  if  the  model  is  substituted  by: 


Tr+V, 

Gr  Gl 

Qr  Qi 

bR  +  b, 

_L 

nR-n, 

i-(rR-r,)\ 

—Gi  Gr 

-Qi  Qr. 

~(bR  -  bj). 

I 

~{nR  +  n;). 

The  derivation  is  as  can  be  simplified  by  noticing  that: 


rR+r, 
~(?R  ~  D) 


rR 

1 

0 

i 

rR 

-D 

-i 

0 

r,. 

and 


r°  i] 

r  A  B 

0 

A 

B 1 

i-i  oJ 

l-B  A- 

-I 

oJ 

—B 

a\ 

Using  these  expressions,  it  is  easy  to  compute: 
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r° 

rR 

Li 

oJ 

fi- 

0  I 

\GR 

G/l 

0 

/ 

0 

/ 

qr 

<2/1 

0 

/ 

0 

/ 

-I  0 

l-G/ 

-/ 

0 

-/ 

0 

l -Q / 

(?r. 

-/ 

0 

-/ 

0 

+ 


0  / 
-/  0 


b. r 

A. 

nR 


0 

I] 

rR 

-I 

oJ 

ft- 

gr  g, 

qr 

Qi 

'  \  ' 

_L 

'  n7  ' 

Gr 

l-Qi 

Qr. 

-  bRm 

1 

■~nR. 

and  therefore, 


rR  +  r, 

Gr  Gi 

qr 

Qi 

bR  +K7 

nR-n, 

-(rR-r;)J 

—Gl  Gr_ 

-Qi 

Qr. 

-ik-bjl 

1 

~(nR  +  n,)_ 

□ 


Approved  for  Public  Release;  Distribution  Unlimited. 
31 


C.  Key  Propositions  and  Proofs 


Proposition  1:  The  product  of  a  diagonal  matrix  DLxL 
commutates  as  follows: 


Dr  —  diag(r )  d 


Proof: 


Di,j  dt  S[j 

C Df)i=  Stj)  Tj 

i 

(D  r)i  —  rj)  dt 

i 

(D  r)i  —  (  diag(f)  d). 


—  diag(d )  and  a  vector  rLxl 


□ 
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Proposition  2:  The  product  of  the  block  matrix 


with  diagonal  matrices  GR  =  diag(gR )  and  G1  =  diag(gj) 


and  a  block  vector  R2lxi  with  entries 


R 


2Lxl 


rR+r, 
i-(rR-r,)  \ 


commutates  as  follows: 


diag(r ',) 
—diag  ( rR ) 


diag(fR)  gR  +  g, 
diag  ft)\  gR  -  g, 


Proof: 


Gr  —Gl 

Tr+T, 

fti  G R  - 

~(fR  +  r,)GR—(fR  —  r,)G, 
_(fR  +  rI)GI-(fR 


The  desired  result  is  obtained  by  applying  Proposition  1: 


diagftR  +  r,)gR-diagftR  -  f{)gi 
.diag ft  +  r^ft-diagift  -  r,)gR. 


diagif \) 
—diag  ( rR ) 


diag  ft)  1  gR  +  g, 
diag  ft)  \  gR  -  gft 
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Proposition  3:  The  average  of  a  function  f(C  ■  b )  over  C  and  b  where  C  G  {±l}Lxt/ 
and  b  G  {+l}t/xl  is  given  by: 

c  b  C 

Proof: 

We  define  zt  j  —  ct  j  bit  then 

c  b  z  b 

-2 

The  number  of  combinations  of  bis2u ,  therefore: 

=  5tfX/(z'T)' 

Z 

□ 
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Proposition  4:  Let  b  —  bR  +  ib1  with  bE{—  1, +1,  — 1,1  }l/xl,  bR  —  Re{b)  and 
bj  —  Im{b}.  Also,  letb+  G  {±l}yxl  and  b_  G  {±1}UX1,  then 

+  b,))f  (F(S„  -  b,))  = 

bR  bi  u  v 


Proof: 


The  set  of  all  possible  combinations  of  bR  and  b,  is  shown  in  the  next  table. 


(5«), 

(*l)f 

(*+),-  =  (*«  +  5/), 

(*-),  =  (*»  -  s,), 

+  1 

0 

+1 

+1 

0 

+1 

+1 

-1 

-1 

0 

-1 

-1 

0 

-1 

-1 

+1 

When  bR  +  bj  —  1,  bR  —  bt  can  be  either  -1  or  +1.  Also,  when  bR  +  b,  —  —1,  bR  —  h/ 
can  be  either  -1  or  +1.  Therefore,  for  summation  purposes,  bR  +  b,  is  statistically 
independent  from  bR  —bh  In  addition,  the  range  of  summation  is  in  {±1}. 

Z f  + s,)) f  ^r(Ss  ■ s,)) = 

bR  b , 


=  ^ZZ/(«IA)/0?TS-) 

b+  b _ 


=  ^j^/(aT£+) 

b+  b _ 

□ 
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Proposition  5:  Letx  and  w  be  vectors  in  RLxl.  The  average  of  the  exponent  of  the 
dot  product  between  x  and  u  over  the  vector  elements  xt  —  ±at  is  given  by: 


\  ^  ^  eUT*  =  J^coshCuf  di). 

Xo=+a0  xi.-i=±ar.-i  i= 0 


xL_1=±aL_ 


Proof:  This  property  can  be  proven  by  induction.  Consider  vectors  x'  and  u'  be 
vectors  inEL+lxl.  The  dot  product  of  vectors  x'  and  u'  can  be  constructed  by 
multiplying  by  the  exponential  term  eUvXL  on  both  sides  of  the  equation: 


\  I  -I  I 


L- 1 


,ul  x+uL-xL  _  „UL 


'Xl  J~|  cosh(Uj  ■  af) 


x0-±a0  xL-1-±aL-1 

Then,  averaging  over  the  variable  xLresults  in: 


1  =  0 


L-l 

\  1  4  1  eux+ULxL  2  I  g  UL-XL  n  COSh(llj  ■  di) 


x0=±a0  xL-±aL 


xL=±aL 


1  =  0 


L-l 

=  cosh(eUL'*L)  J- 1  cosh(iq  ■  cq) 
i= o 
L 

=  J” |  cosh(iq  ■  di) 

i= 0 

Therefore,  the  formula  applies  to  vectors  in  Ri+1. 


□ 
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Proposition  6:  The  TSC  is  invariant  to  permutations  and  sign-inversions  of  rows 
and  columns. 


Proof:  A  permutation  matrix  T  has  the  following  property: 

-?l  =  I  and  P2  ■  Pj  =  / 

Using  our  definition  of  a  TSC  for  a  binary  antipodal  matrix, 

t2(C)  -  \\ChC\\2f  -L-U, 

we  have 

t2(Pi  CP2)  =  ||P2TCrP1TP1CP2|||.  —  L-U 
=  ||CrP1rP1CP2P27’|||  -  L-  U 
=  \\CtC\\2f-L-U 
=  T  2(C). 

□ 
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The  following  definitions  are  related  to  Proposition  7. 


r  -.Lxn 

Definition:  The  set  S 3  is  the  subset  of  real  antipodal  matrices  C  —  (q j  j  with: 

u-i 

ai  =  ^  Cij  >  0. 

7=0 

r  -,1x1/ 

Definition:  S3 1  is  the  subset  of  real  antipodal  matrices  C  —  {q;- j  G  S3  with: 

t2(C)  =  t. 

Definition:  S^+  ^_  t  is  the  subset  of  pair  of  matrices  (Z+,Z_)  with  Z+  =  {±l}Lxt/ 
and  Z_  =  {±l}ixU  with: 

u-i 

(a+)i  =  >  0 

7=0 

£7  —  1 

(a-)i  =  ^(Z-)ij  >  0 
7=0 

2  / (z+  +  z_)  _  .(Z+-Z_)\ 

1  = 
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Proposition  7.  The  limit  of  aLUgQ3)as  /?  approaches  infinity  depends  on  the  feature 
vectors  as  follows: 


lim  aLiUiS  (/?) 

p  — >oo 


0  for  a  *  a* 


2Sum(ai=0') 


S fit  ■ 

u’Lmin 


Xa 


2l 


>  0  for  a  —  a* 


2Sum{a^— 0)  |  ^Xmin  I 


Proof:  The  expression  can  be  converted  in  a  polynomial  rational  function  by 
substituting  e~P  ->  y  and  divide  by  the  exponent  with  the  minimum  Total  Square 
Correlation  yTmin  in  both  side  of  the  equation. 


aL,£/,a(/?) 


2Count(a 


£— o)  y 


2l 


Sa  2^oun^iai=0)  ZtZzes^y1 


2Count(a.i=0 ) 


=o)  ^  t  y 


Sa: 


2l 


:  2 Count (a^= 0 J  2 1  2jZE  ,y;  f  T 

t  =  t(Z)  >  tmin  >  0 
a  =  abs(Z  ■  1) 

Since  t  —  tmin  >  0,  taking  the  limit: 


t-t-n 


lim  yt  tmin  — 
y->o  fl 


0  fort*  tmin 
for  t  =  tmin 


eliminates  all  terms  with  Total  Square  Correlation  terms  greater  thanTmjn.  The 
result  is: 


aL,U,a(P ) 


z  y 

2Count{ai= 0)  ZaZE  Sa,tmi 

XX 


'  2^ount(aj=0)  ZuZE  S-^  t 


>  o 


The  summation  over  all  the  terms  in  the  set  is  just  the  cardinality  of  the  set, 
therefore: 


aL,U,a(P}  ~ 


2Cownt(aj-0)  |  Q->tmin  I 
,  2L  |  „ 


■>a  2Count{aj=0)  \  ^ 0.,tmin\ 
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□ 


D.  Orthogonal  Relationships 

Let  G  —  GR  +  G,  be  a  diagonal  matrix  with  diagonal  elements  given  by 

G  —  diag(gR  +  i  g ,) 

and  gR  +  n  g,  G  {+1,  ±i}Lxl. 

Then  the  following  equations  are  true: 

GrG,  —  diag(gR )  diag(g ,)  —  0 
GrGr  +  Gj  Gi  —  I. 


Proof: 

Whenever  g  takes  a  real  value,  the  imaginary  part  is  zero  by  definition. 

(Hr  E  ^ 

Im{(gR  +ig1)i)  =  0 

(9i)t  =  0 

Whenever  g  takes  a  imaginary  value,  the  real  part  is  zero  by  definition. 

(, gR  +  fl  g[)i  E  Imaginary 
Re{(gR  +5  #/)*}  =  0 
(Hr)i  =  0 

Therefore,  the  product  GRG,  is  a  zero  matrix. 

GrGi  —  diag(gR )  •  diag(gj') 

=  diag({(jjR)i(gI')i}LX1') 

—  diag({  0  }ixl) 

=  0. 
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By  the  definition  of  gR  and  gh  we  have  that: 


9r  +  9i  =  ±1- 

Using  the  previous  results,  we  have  that  the  identity  is  formed  by  GRGR  +  U/ G,. 

I  =  diag(gR  +  g,)  ■  diag(gR  +  g,) 

I  =  diag(gR  +  g,)  ■  diag(gR  +  g,) 

I  —  diag(gR)T  ■  diag(gR )  +  diag(gRY  ■  diag^g,) 

+diag(g,Y  ■  diag(gR )  +  diag(gjY  ■  diag(gj) 

I  —  diag(gRY  •  diag(gR )  +  diag(g,Y  •  diag(g ,) 

I  ~  GrGr  +  Gj Gj 
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8.  Acronyms 


ALRT 

CDMA 

BPSK 

QPSK 

LRT 

MAP 

TSC 

t2() 

y 

E 

A(f  |  JO 
A(r \‘K,  C,  b) 
C 

pLxU  f 
L  >  lLxU 

C*,i 

c  r 

Cr 
Cl 
c * 

b 
b r 
bi 

C 9)i 

aL,U,a 

a * 
a* 

Y/(a) 


Average  Likelihood  Ratio  Test 
Code  Division  Multiple  Access 
Binary  Phase  Shift  Keying 
Quadrature  Phase  Shift  Keying 
Likelihood  Ratio  Test 
Maximum  A  Posteriori  Probability 
Total  Squared  Correlation 
Total  Squared  Correlation  function 
Chip  Signal-to-Noise  Ratio 
Energy  per  symbol  per  user 
Likelihood  function 
Conditional  likelihood  function 

Spreading  matrix 

Matrix  with  specified  rows  (L)  and  columns  ([/) 

Column  vector  i  of  matrix  C 

Hermitian  operator,  i.e.,  the  complex  conjugate  transpose 

Real  part  of  the  matrix  C 

Imaginary  part  of  the  matrix  C 

TSC  optimal  spreading  matrix 

Information  vector 

Real  part  of  the  vector  b 
Imaginary  part  of  the  vector  b 
Element  i  of  vector  g 

Coefficients  of  the  CDMA  likelihood  function 
Non-vanishing  alpha  coefficients 
Feature  vector  for  CDMA  detection 
Summation  over  all  possible  vector  coefficients 
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I  AO 

Summation  over  all  possible  matrix  coefficients 

c 

I 

Identity  matrix 

0 

Zero  matrix 

1 

Row  Column  vector  with  elements  equal  to  1 

diag  ( g ) 

Diagonal  matrix  with  elements  (g)  t  in  the  diagonal. 

E{  } 

Expectation  operator 

R 

Set  of  real  numbers 

\s\ 

Cardinality  of  a  set  S,  i.e.,  the  total  number  of  elements  in  the 
set 

{±1} 

Set  of  antipodal  symbols 

{±1,  ±a} 

Set  of  QPSK  symbols 

Approved  for  Public  Release;  Distribution  Unlimited. 
43 


